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Abstract 

With  a  lot  of  experiments  on  explosion,  the  dynamic  stability  of  grotto  and  sprayed  anchor 
strutted  grotto  were  studied  under  effects  of  seismic  wave  of  explosion,  and  a  formula  was  put 
out  for  calculating  the  stability.  The  results  of  calculation  fitted  results  of  testing  well. 

1.  Introduction 

It  is  usually  needed  to  evaluated  the  stability  of  tunnels  and  galleries  in  structures  of  mines, 
railroads  and  hydro-electric  engineerings  under  effects  of  dynamic  loading,  and  the  safety 
distance  needed  to  be  determined.  The  dynamic  loads  result  mainly  from  explosion  work  and 
accidents.  When  the  stability  of  tunnels  and  galleries  are  evaluated  under  effects  of  seismic 
wave  of  explosion,  it  should  be  considered  the  dynamic  loads  from  explosion  as  well  as  static 
loads  from  rocky  soil.  The  quantities  of  dynamic  loads  are  related  to  amplitude  and  lasting 
period  of  seismic  wave  of  explosion  propagated  among  ground. 

2.  The  dynamic  strength  of  rocks 

It  is  well  known  that  under  effects  of  dynamic  loads,  the  limit  range  of  rock  strength  increases 
with  the  increasing  of  loading  rate.  The  increased  value  of  strength  relates  to  the  nature  of 
rocks  and  loading  rate.  The  values  of  granite  and  marble  are  listed  in  table  1. 
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Table  1.  The  variations  of  compressive  strength  of  rocks  with  loading  rate 


Rock  Loading  rate  Com.  strength  Loading  rate  Com.  strength  c^/ci 
Vi  kg/crtWsec  Ci  kg/cm2  V2kg/cm2/sec  rr2  kg/cm2 


Granite  5  1220  20  2000 

Marble  5  500  30000  980 


1.64 

1.98 


Table  1.  The  variations  of  compressive  strength  of  rocks  with  loading  rate 


All  of  compressive  strength  of  rocks  increase  with  increasing  of  loading  rate,  but  increased 
value  is  different  for  each  rock.  The  compressive  and  flexural  strength  and  elasticity  modules 
of  rock  increase  with  a  logarithmic  function  of  loading  rate. 

3.  The  dynamic  strength  of  rock  mass 

Because  cracks  and  cleavages  exist,  the  strength  of  rock  mass  is  lower  than  that  of  rock.  The 
discount  factor  is  about  0.80-0.90.  When  the  loading  rate  is  102~104kg/cm2/sec,  the 
compressive  strength  of  rock  increases  by  1. 1  6-1  .43  times  and  flexural  strength  increases  by 
1  .24-1  .48  times.  Because  the  damage  of  rock  mass  of  tunnel  is  controlled  by  tensile  strength, 
the  dynamic  strength  of  rock  is  taken  as  1.30-1.40  times  of  static  strength.  So,  the  formula  for 
calculating  the  dynamic  strength  of  rock  mass  is: 

gd  =  Kdgp  (1) 

Where:  gd  is  dynamic  tensile  strength  of  rock  mass  (kg/cm2);  GP  is  ultimate  static  tensile 
strength  of  rock  (kg/cm2);  KD  is  increasing  factor  of  dynamic  strength  of  rock  mass. 

If  the  surface  rocks  of  grotto  are  stable,  the  surface  are  usually  sprayed  with  Scm  thickness 
concrete  and  K  is  among  1.04-1.26.  If  the  surface  rocks  of  grotto  are  unstable,  the  surface  are 
usually  strutted  by  rollbolts  and  then  sprayed  with  Scm  thickness  concrete  and  K  is  among 
1.30-1.40. 


4.  The  calculation  of  stability  of  tunnels  under  effect  of  seismic  wave  of  explosion 

4.1.  Concentrated  factors  of  dynamic  and  static  stress  of  tunnel  without  lining 

Under  effect  of  plane  wave  of  explosion,  the  concentrated  factors  of  dynamic  and  static  stress 
of  tunnel  without  lining  are  listed  in  table  2. 


Table  2.  Concentrated  factors  of  dynamic  and  static  stress  of  tunnel  without  lining* 


Pattern  of  Position  Concentrated  factor  of  dynamic  stress  Concentrated  factor 
tunnel  - - -  of  static  stress 

Analytical  D/X=0. 12~0. 15  Numerical  D/%=0.50 

Straight  wall  Arch  3.25~2.25  3.00  3.25 

and  round  arch  Wall  2.00-1.65  1.80  1.50 


*D  is  diameter  of  tunnel;  X  is  wave  length 


Table  2. 

Concentrated  factors  of  dynamic  and  static  stress  of  tunnel  without  lining* 


4.2.  The  calculation  of  stability  of  tunnel  under  effects  of  seismic  wave  of  explosion 

Under  effects  of  seismic  wave  of  explosion,  the  stable  condition  of  rock  mass  of  tunnel  without 
lining  is  that  the  sum  of  static  stress  of  mountain  and  dynamic  stress  of  seismic  wave  of  explosion 
is  smaller  than  or  equals  to  the  dynamic  strength  of  rock  mass;  that  is: 

o  =  oCT  +0DT  <[oD]  (2) 

Where:  o  is  the  total  stress  among  rock  mass  (kg/cm2);  [aD]  is  allowable  dynamic  strength  of 
dynamic  strength  of  rock  mass  calculated  from  equation  (1)  (kg/cm2);  oDX  is  dynamic  stress 
of  rock  mass  under  effect  of  seismic  wave  of  explosion,  (kg/cm2). 

oDT  =  (l/K0)(KGy/g)CeV  x  103  (3) 

Where:  y  is  unit  weight  of  rock  (kg/cm2);  V  is  vibrating  velocity  of  rock  particles  in  periphery 
of  tunnel  under  effects  of  seismic  wave  of  explosion  (cm/sec);  Ce  is  elastic  velocity  of 
longitudinal  wave;  g  is  acceleration  of  gravity  g=9.81  m/sec2;  K0is  reflection  factor.  Tested 
results  showed  that  if  the  tunnel  is  acted  by  incident  seismic  wave  of  explosion,  the  factor  is 
K0=2;  if  the  tunnel  is  acted  by  reflected  seismic  wave  of  explosion,  the  factor  is  K0=  2. 

4.3.  The  calculation  of  critical  vibrating  velocity  of  rock  particles 
From  equation  (2),  the  critical  vibrating  velocity  of  rock  particle  is: 

Ve  =  K0(Kdgp  -  oCT)/(KGyCe)  x  1  0"3  (4) 

Where:  Ve  is  critical  vibrating  velocity  of  rock  particle.  When  the  critical  vibrating  velocity  of 
rock  particle  is  calculated  to  rock  mass  in  elasticity  zone,  the  elastic  velocity  of  longitudinal 
wave  Ce  is  taken;  when  the  critical  vibrating  velocity  of  rock  particle  is  calculated  to  rock 
mass  with  cracks,  the  elasto-plastic  velocity  of  longitudinal  wave  Cp  is  taken,  and  it  is  taken 
Cp=Ce/2  if  there  is  no  tested  data  of  Cp. 

4.4.  The  calculation  of  critical  vibrating  velocity  of  rock  particle  for  collapsed  rock  mass 

The  tested  results  showed  that  when  the  properties  of  tunnel  structure  change  to  plasticity, 
cracks  appear  in  arch  and  in  boundaries  between  arch  bottom  and  side  wall.  With  the 
continuous  effect  of  explosion  wave,  the  deformations  of  arch  and  side  wall  increase,  but  the 
stress  in  rocks  doesn't  increase  and  the  tunnel  appears  a  unloading  effect.  Under  effect  of  a 
normal  big  explosion,  the  self-vibrating  frequency  of  rock  is  10-15  Hz  and  the  seismic  wave 
of  explosion  lasts  0.4-0. 6  sec.  If  the  tunnel  is  taken  as  a  structure,  the  unloading  factor  of 
tunnel  in  plasticity  zone  is  solved  from  theory  calculation,  thus  the  critical  vibrating  velocity 
of  rock  particles  in  collapsed  rock  mass  is: 


VP  =  K0(KDoP  +  oCT)g/(KGyCP)  •  (1/KJxlO'3  (5) 


With  tested  data  of  deformation  of  tunnel  caused  by  effect  of  seismic  wave  of  explosion,  the 
unloading  factor  Kz  is  solved  theoretically.  If  unloading  factories  Kz=0. 80-0.65,  rock  mass 
appears  partial  collapsing  and  the  collapsed  volume  is  usually  smaller  than  Im3.  If  unloading 
factor  is  Kz  =  0.50-0.35,  the  rock  mass  appears  large  area  collapsing. 

4.5.  Comparison  between  results  of  calculation  and  practical  tests 

Known  a  round  arch  and  straight  wall  tunnel  without  lining,  rollbolt  and  sprayed  concrete. 
The  span  is  L=3m,  height  is  H=3m.  The  rock  is  granite  with  macro  crystalline  (weathering). 
The  properties  of  rock  are  tested  as:  Pope's  factor  f=4~6;  unit  weight  y -2. 64t/m  ;  elastic 
velocity  of  longitudinal  wave  Ce=2060m/sec;  dynamic  elastic  modules  E=0.928xl05kg/cm2; 
Poisson  ratio  (I  =  0.30;  internal  friction  angle  (])  =  41°06';  static  ultimate  tensile  strength  ap  = 
23.60kg/cm2. 

Chosen  f=5,  Ce=2060m/sec,  Cp=1030ni/sec,  Kp=1.15,  K0=2,  Kz=0.65  for  partial  collapsing 
and  Kz=0.35  for  large  area  collapsinng,  the  results  of  calculation  is  listed  in  table  3. 


Table  3.  A  comparison  of  critical  vibrating  velocity  of  rock  particles 
from  calculation  and  from  practical  tests  cm/sec 


No  damage 
LV1] 

Slight  cracking 
[V2] 

Partial  collapsing 
[V3] 

Large  area  collapsing 
[V4] 

Tested 

30 

30-50 

50-100 

100-200 

Calculated 

30.36 

30.36-60.73 

60.73-93.42 

93.42-173.48 

Table  3.  A  comparison  of  critical  vibrating  velocity  of  rock  particles  from 
calculation  and  from  practical  tests  cm/sec 


The  calculated  results  fit  the  tested  results  well  and  the  maxiinum  error  is  smaller  than  20%. 

5.  Applications  in  engineering 

5.1.  Land  form 

If  the  thickness  of  rock  mass  at  minimum  resistant  line  of  grotto  is  smaller  than  50  times  of 
explosive  diameter,  the  rock  mass  is  called  gentle  slop;  if  the  thickness  of  rock  mass  is  bigger 
than  50  times  of  explosive  diameter,  the  rock  mass  is  called  steep  slop. 

When  a  explosive  in  grotto  explosion,  a  projection  of  rock  over  the  grotto  takes  place  if  the 


grotto  is  in  a  gentle  slop  and  projection  doesn't  take  place  if  the  grotto  in  a  steep  slop. 

5.2.  Geology 

The  classification  of  rock  is  made  according  to  features  of  structure  and  is  listed  in  table  4. 


Table  4.  Classification  of  structure  of  rock  mass 


Classification  Structure  features  Comp,  strength  Elastic  velocity  n 

of  rock  of  longitudinal 

kg/cm  wave  m/sec 


Concordant 

structure 


Massive 

structure 


Fragment 

structure 


Rock  mass  is  a  whole  or  a 
giant  layer,  extramly  undeveloped 
joints,  no  dominate  strructure 
plane  Bo*=l— 2,  M<0.5 
Rock  mass  is  a  massive  or  thick 
layer  structure,  undeveloped  joints 
most  of  structure  planes  are  joint 
plane  and  closed  (such  as  psephyte) 
Bo=2~3,  M=0.5~2 
Rock  mass  is  a  less  thick  layer  or 
massive  structure,  developed  joint, 
most  of  structure  planes  are  jiont 
planes.  Bo=3— 4.  M=2~5 


>300 


>200 


>100 


>4000  >0.85 


3000-4000  0.85-0.6 


2000-3500  0.6— 0.3 


*Bo  is  joint  data;  M  is  number  of  joints  in  one  meter;  n=(Cv/Ce)2;  Cv  is  longitudinal  wave  velocity  of 
rock  mass  m/sec.  Ce  is  longitudinal  wave  velocity  of  rock  m/sec. 


Table  4.  Classification  of  structure  of  rock  mass 


5.3.  The  critical  vibrating  velocity  of  wall  rock  with  different  degree  of  damage 


The  velocity  are  listed  in  table  5 . 


Table  5.  Critical  vibrating  velocity  of  wall  rock  of  grotto  [VI],  [V2],  [V3],  [V4] 


Rocks  Unit  weight  Comp,  strength  Tens,  strength  No  damage  Slight  damage  Intermediate  damage  Serious  damage 


— 

(t/m3) 

(kg/cm^) 

(kg/cm^) 

(cm/sec) 

(cm/sec) 

(cm/sec) 

(cm/sec) 

Hard 

2.60-2.70 

750-1 100 

21-34 

27 

54 

82 

153 

rock 

2.70-2.90 

1100-1800 

34-51 

31 

62 

96 

178 

2.70-2.90 

1800-2000 

51-57 

36 

72 

111 

209 

Soft 

2.00-2.50 

400-1000 

11-31 

29 

58 

90 

167 

rock 

2.00-2.50 

1000-1600 

34-45 

35 

70 

107 

199 

Note:  1.  The  data  in  this  table  are  applicable  for  grottoes  in  gentle  slop;  if  grottoes  in  steep  slop, 
the  data  needed  to  be  multiplied  by  2. 

2.  If  the  hole  of  explosives  is  parallel  or  oblique  or  perpendicular  to  a  grotto,  the  critical 
vibrating  velocity  of  wall  rock  [VI],  [V2],  [V3]  and  [V4]  need  to  be  multiplied  by  1.0,  1.2 
and  1.4,  respectively. 

3.  The  data  in  the  table  are  applicable  for  rocks  with  concordant  structure  and  the  data  need 
to  be  multiplied  by  0.9  or  0.8  for  rocks  with  massive  structure  and  fragment  structure, 
respectively. 

Table  5.  Critical  vibrating  velocity  of  wall  rock  of  grotto  [VI],  [V2],  [V3], 

[V4] 


5.4.  The  determination  of  support  pattern  of  wall  rock 

According  to  the  relation  between  perpendicular  vibrating  velocity  of  seismic  wave  of 
explosion  and  critical  vibrating  velocities  of  wall  rock  of  no  damage[Vl],  slight  damage[V2], 
intermediate  damage[V3]  and  serious  damage[V4],  the  support  pattern  of  wall  rock  is 
determined: 

A.  if  Vv<[Vl],  sprayed  with  #200  plain  concrete  of  thickness  50  cm; 

B.  if  [VI]<Vv<[V2],  sprayed  with  #200  plain  concrete  of  thickness  80  cm; 

C.  if  [V2]<Vv<[V3],  net  and  concrete;  a  250x250  mm  net  made  of  8mm  steel  is  placed  on  the 
surface  of  wall  rock  and  then  sprayed  with  #200  plain  concrete  of  thickness  80~  100mm; 


D.  if  [V3]<Vv<[V4],  net,  concrete  and  roltbolt;  a  250x250  mm  net  is  placed  on  surface  and 
sprayed  with  -200  plain  concrete  of  thickness  80-100  mm,  and  -16x2000  mm  rollbolts  are 
installed.  The  rollbolts  are  made  of  mortar  and  arranged  with  a  distance  of  2000x2000  mm  to 
each  other. 


6.  The  calculation  of  perpendicular  vibrating  velocity  of  seismic  wave  of  explosion 

Based  on  analysis  of  tested  data  from  explosion  on  ground,  mine  explosion  including  standing 
shot,  long  hole  volley  fning,  long  hole  short-delay  blasting,  directional  explosion  and  internal 
explosion  of  tunnel,  empirical  formulae  for  calculating  the  perpendicular  vibrating  velocity  of 
rocks  in  various  geological  conditions  are  listed  in  table  6  and  drawn  in  fig.  1  and  2. 


Table  6.  Empirical  formulae  of  perpendicular  vibrating  velocity  of  rock 
particles  under  effects  of  seismic  wave  of  explosion 


No  Pattern  of 

Explosion  conditions  and 

Grological 

Vv=k(Ql/3/R)a 

explosion 

explosive  quantity 

condition 

k 

a 

1  ground 

Central  charged 

Granite 

98.76 

1.37 

explosion 

Q=l,  3,  5,  10,  15,  14,1  OOt 

2  Unshelted 

A.  once  delayed  blasting 

Diabase 

804 

2.24 

big  explosion  Q=9320t 

Disbase 

630 

2.80 

©Standing 

B.  Vollery  blasting  Q=1000t  Diabase 

206 

1.81 

shot 

Q=534t  Mctamorphic  rock 

180 

1.47 

Q=  lll~178t  Metamorphic  rock 

79 

1.39 

Phyllite 

82.5 

1.32 

Q=20t  Mica-quartz  schist 

152 

1.56 

Phyllite 

156 

1.93 

Q=305t 

Diabase 

718 

2.40 

Granite  and  marble 

150 

2.00 

©Long-hole 

Q=200t 

Marble 

77.6 

2.33 

vollery 

Q=103t 

Quartz 

624 

2.41 

blasting 

Q=8~l  4t 

Limestone 

125.7 

1.63 

Limestone 

130 

1.80 

Limestone 

140 

1.80 

Limestone 

200 

1.80 

Limestone 

340 

1.80 

©Long-hole 

A  .  6  pieces  Q=45.9 

Gneiss 

180 

1.83 

shoet  ddelay 

Gneiss 

116.2 

1.73 

blasting 

B.  10  pieces  Q=4.23t 

Marble 

378 

1.60 

Marble 

107 

1.50 

C.  1 0  pieces  Q=4.74t 

Primary 

130 

1.70 

D.  10  pieces 

Quartz 

142 

1.61 

Quartz 

153 

1.60 

©Directional  A.  Total  Q=13394t 

Sandstone 

240 

2.00 

blasting 

B.  Total  Q=503t 

Diabase 

115 

2.00 

3  Internal 

Linear  charge  Q=8~15t 

Granite 

99.6 

1.72 

explosion  of 

Granite 

111.2 

1.92 

tunnel 

Granite 

591.4 

2.30 

Granite 

90.8 

1.82 

Note:  Vv  is  perpendicular  vibratting  velocity  of  seismic  wave  of  explosion  (cm/sec); 
Q  is  total  quantity  of  explosion;  R  is  distance  between  centre  of  explosion  and 
testing  point  (m). 


Table  6.  Empirical  formulae  of  perpendicular  vibrating  velocity  of  rock 
particles  under  effects  of  seismic  wave  of  explosion 


Fig.  1. 

Distribution  of  perpendicular  vibrating  velocity  of  seismic  wave  of 

explosion 


Scale  quantity  of  explosive  (kg1/3/m) 


Fig.l.  Distribution  of  perpendicular  vibrating  velocity  of  seismic  wave  of  explosion 


Fig.2.  Distribution  of  velocity  of  seismic  wave  of  internal  explosion  of 

tunnel 


0.1  0.2  0.3  0.6  1  2  3  4  5 

Scale  quantity  of  explosive  (kg1/3/m) 

Fig.2.  Distribution  of  velocity  of  seismic  wave  of  internal  explosion  of  tunnel 


Fig.  1  and  fig.2  show  relations  between  perpendicular  velocity  of  rocks  and  scale  quantity  of 
explosiveQ1/3/R  under  various  conditions  of  explosion. 


7.  Safety  distance  of  tunnel  without  lining  under  effect  of  seismic  wave  of  explosion 

With  the  empirical  formulae  and  tested  data  in  table  6,  the  calculating  formula  of  safety 
distance  is  derived: 

R  =  l/([V]/k)(1/ct)  Q1/3  (6) 

Where:  R  is  safety  distance  (m)  of  tunnel  without  lining  under  effects  of  seismic  wave  of 
explosion;  a,  K  are  factors  determined  from  tested  data  in  table  6  or  in  fig.  1~2.;  Q  is 
calculated  explosive  quantities:  total  quantities  of  standing  shot,  the  maximum  quantity 
among  each  delay  explosion  and  the  maximum  quantity  in  each  period  of  short-delay 
blasting;  [V]  is  critical  vibrating  velocity  of  rock  particles  (cm/sec)  and  calculated  with 
equation  (4),  (5)  or  with  formulae  in  table  5. 


8.  Conclusions 


1.  On  basis  of  balance  between  dynamic  strength  of  rock  and  the  sum  of  dynamic  and 
static  stress  acted  on  tunnel,  the  formulae  for  calculating  critical  vibrating  velocity  of  rocks 
are  derived  when  tunnel  without  lining  un  elastic  state  appears  cracking,  partial  collapsing  and 
large  area  collapsing  under  effect  of  seismic  wave  of  explosion.  The  results  of  calculating  fit 
the  tested  results  well. 

2.  With  the  critical  vibrating  velocity,  the  support  pattern  of  rollbolt  and  sprayed  concrete, 
and  safety  distance  of  stability  of  grotto  are  determined. 


